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I. INTRODUCTION 



A precise measurement of the W boson mass {My/) is of fundamental importance in 
testing the standard model (SM) and in constraining new physics. Future experiments 
at LEP2 and the Tevatron aim at a precision of about 40 MeV and 20 MeV for 
My/, respectively. In order to be able to achieve a precision of O(10 MeV) for My/ in 
hadron collider experiments, it will be necessary to control not only the higher order QCD 
corrections, but also the electroweak radiative corrections. 

At 0(a n ), W decay with collinear emission of photons from a final state charged lepton 
of mass me gives rise to terms which are proportional to (a/ir) n \og n (M v/ /mj), in n-photon 
exclusive rates. These terms can have significant effects in determining My/- At 0(a), 
photon emission shifts the W mass by [|J 

AM w ^-^-T w , (1) 

where Y w is the total width of the W boson, and 

(2) 

When finite detector resolution effects are included, IV is replaced by the uncorrected 
experimental width. For CDF, the resulting shift is —168 MeV (—65 MeV) in the muon 
(electron) channel [f|. Similar results are found in the D0 W mass analysis ||. In the 
electron case, understanding the exact amount of energy lost in photon bremsstrahlung is 
also important to determine the E/p distribution, which is used by CDF to determine the 
energy scale of the central electromagnetic calorimeter ||. When aiming for a precision 
of O(10 MeV) in My/, it is therefore important to control the effects of multiple photon 
radiation in W events. 

Here, we present a calculation of two photon radiation in the W —>■ \iv channel at the tree 
level, including the finite mass of the muon which regulates the collinear singularity originat- 
ing from photons radiated off the final state muon. Initial and final state bremsstrahlung, 
together with finite W width effects are taken into account in our calculation. 



II. CALCULATIONAL DETAILS 



The cross section for pp — ► £"^77 has been calculated in the limit of massless fermions 
in Ref. M and [EL using the helicity amplitude technique of Ref. J5|. To evaluate the matrix 
elements for a massive final state charged lepton, we use the MADGRAPH [nj package, 
which automatically generates the SM matrix elements in HELAS format fil l. 

However, in order to maintain electromagnetic gauge invariance in presence of finite 
W width effects, the W propagator and the WW'j and WW'yy vertex functions in the 
amplitudes generated by MADGRAPH have to be modified f8j,[12||. Finite width effects are 
included by resumming the imaginary part of the W vacuum polarization, Ily^{q 2 ). The 
transverse part of Uw{q 2 ) receives an imaginary contribution 



Imn^(g 2 ) = g 



2 rv 



M 



w 



(3) 



while the imaginary part of the longitudinal piece vanishes. The W propagator is thus given 
by 



q 2 - Mfr + iq 2 j w 
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(4) 



with 
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(5) 



A gauge invariant expression for the amplitude is then obtained by attaching the final state 
photons to all charged particle propagators, including those in the fermion loops which 
contribute to Uw(q 2 )- As a result, the triple and quartic gauge vertices and ) 

are modified to 



(6) 
(7) 



Due to the mass singular terms which arise from the fi'j collinear region, it is nontrivial 
to numerically compute the ^^77 cross section. To optimize the Monte Carlo integration, 



we choose the logarithms of the invariant masses of the /r/i, the ^72, and the ^7172 system 
as integration variables: 

dPS± ~ ds rflogm(/i7 1 7 2 ) dlogm(/i7i) dlogm(/x7 2 ). (8) 

The resulting cross section is checked for gauge invariance, and is numerically stable 
when the photons are collinear with the muon. 

For pp — > ez/77, due to the smaller electron mass, the singularities for both photons and 
the electron being collinear, are much more severe than in the muon case. The method used 
in pp — > /iz/77 to perform the phase space integration is not sufficient to obtain a stable 
numerical answer for pp —>■ ewyy. 



III. RESULTS AND DISCUSSION 



For our numerical simulations we use the MRSA set of parton distributions |L3| , M\y = 
80.22 GeV, and a = 1/128. The factorization scale is fixed to the parton center of mass 
energy, y§. To simulate detector response, we impose the following acceptance cuts on 
transverse momenta and pseudorapidities: 

PtQj.) > 25 GeV, i> T > 25 GeV, (9) 

\ V (H)\ < 1.0, M 7 )|<3.6. (10) 

If the two photons have a separation smaller than a critical value Ai?( 7 , 7) = R c , they 
cannot be discriminated in the detector, and thus will be treated as one photon. Taking the 



typical cell size of an electromagnetic calorimeter as a criterion |14}| , we choose R c = 0.14 in 
the following. The results of our simulations are listed in Table |. 

When no separation requirement is imposed on the muon and the photon, approximately 
13% (5.7%) of all W — > [iv events contain one photon with a minimum photon transverse 
energy of £$,(7) = 0.1 GeV (£§.(7) = 1 GeV) @ at 0(a). About 0.8% (0.13%) of all 
events contain two isolated photons with transverse energies larger than £^( 7 ) = 0.1 GeV 



(Et(j) = 1 GeV). If photons are required to be isolated from the muon (AR(fi, 7) > 0.7), 
about 3.8% of the W — > \w events contain one photon with E T {^) > 0.1 GeV, but only 
0.07% have two photons for the same £t(t) threshold (see Table ^p). For larger photon Et 
thresholds, the fraction of two photon events drops very rapidly. 

Although the two-photon tree level calculation presented here is not a full next-to-next 
to leading order calculation of production, we can still gain some important information 
on multiphoton radiation. Knowing the cross section for the lowest order process, <7o, and 
for O(a) one photon emission, a\, one can estimate the probability, P{n), for radiating n 
photons. Neglecting the dependence of the soft photon emission rate on the (hard) quark 
and lepton configuration, the probability is given by [15| 



P(n) = ^^e-<">, (11) 



where < n >= <j\/<j$. Comparison with Table | shows that Eq. ( |TTD provides a useful 
guideline for multiple photon radiation in W events if the minimum photon E? threshold is 
smaller than about 3 GeV. 

In summary, we have calculated two photon bremsstrahlung in W — > [iv events at hadron 
colliders. Our calculation takes both initial and final state photon radiation, and finite W 
width effects into account. We also incorporate finite muon mass effects, which makes our 
calculation valid over the entire phase space region. Approximately 0.8% of all W — > fiu 
events are found to contain two isolated photons with transverse energy £^(7) > 0.1 GeV, 
which roughly corresponds to the tower threshold of the electromagnetic calorimeter of CDF 



and D0 |L4j . Multiple photon bremsstrahlung thus is expected to have a non-negligible effect 
on the W mass extracted from experiment. Unfortunately, our calculation at present does 
not produce stable numerical results in the W — > ev case, due to the more severe nature of 
the collinear singularities in the electron case. It, therefore, should only be viewed as a very 
first step towards a more complete understanding of multiple photon radiation in W boson 
events. 
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TABLES 

TABLE I. Fraction of W — ► jiv events (in percent) containing one or two photons at the 
Tevatron (y/s = 1.8 TeV), (a) with no isolation cut between the photons and the muon imposed, 
(b) requiring Ai2(/z, 7 ) > 0.7. In the two photon case, the two photons are required to fulfill the 
condition AR(j,j) > 0.14. Fractions are obtained by normalization with respect to ctq = 0.346 



nb. 



(a) no isolation cut 




W — > fiwy 


W — ► /xz/77 


0.1 


13.4 


0.77 


0.3 


9.60 


0.39 


1.0 


5.67 


0.13 


3.0 


2.62 


0.024 


10.0 


0.32 


7.6 x 10~ 4 


30.0 


0.013 


1.8 x 10~ 5 


(b)Ai2( M)7 ) >0.7 


^(7) 


W — > (ivy 


W — > /xz/77 


0.1 


3.76 


0.071 


0.3 


2.60 


0.037 


1.0 


1.46 


0.012 


3.0 


0.65 


2.9 x 10~ 3 


10.0 


0.13 


2.0 x 10~ 4 


30.0 


0.012 


7.4 x 10~ 6 
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